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DESIGN, PERFORMANCE, AND CALCULATED ERROR OF A 
FARADAY CUP FOR ABSOLUTE BEAM CURRENT 
MEASUREMENTS OF 600-MeV PROTONS 

Sherwin M. Beck 
Langley Research Center 

SUMMARY 

A Faraday cup was designed and constructed for absolute beam current measure- 
ments of 600-MeV protons produced in the NASA Space Radiation Effects Laboratory. 

The movable Faraday cup assembly has an estimated weight of 616.9 kg. The cup alone 
has an outside diameter of 20.32 cm, an inside diameter of 15.24 cm, a base thickness of 
30.48 cm, a total length of 106.7 cm, and a weight of 235.9 kg. It is constructed of natural 
lead encased in a copper jacket. 

Several possible errors associated with incident proton interactions with various 
parts of the assembly are evaluated. A calculated systematic error of -0.83 percent is 
obtained as the decrease from the true current value. 

From data obtained in calibrating helium-filled ion chambers with the Faraday cup, 
the mean energy required to produce one ion pair in helium is found to be 30.76 ±0.95 eV 
for nominal 600-MeV protons. This value agrees well, within experimental error, with 
reported values of 29.9 eV and 30.2 eV. 

INTRODUCTION 

As part of an experimental program to obtain differential cross sections for the pro- 
duction of charged secondary particles from elemental materials under 600-MeV-proton 
irradiation, a Faraday cup was designed, constructed, and used as an absolute proton beam 
monitor at the NASA Space Radiation Effects Laboratory (SREL). For maximum utiliza- 
tion in the three main experimental areas of the laboratory, the cup was designed as a 
mobile self-contained instrument complete with vacuum pumping equipment, vacuum read- 
out gage, electromagnet, and magnet power supply. 

The design of the cup was based on available information in the literature (refs. 1 
to 6), experience of SREL staff members, and constraints imposed by requiring a mobile 
system. The design was also influenced by the availability of materials and the need to 
complete the cup in a relatively short period of time. In order to produce a mobile cup 



with a 100-percent charge collection efficiency, the arbitrary decision was made to design 
the cup for incident proton beams with circular cross section of 10 cm or less and to make 
the cup base 15 to 20 percent larger than the range of 600-MeV protons. 

This report describes the features of the Faraday cup system and its use as an 
absolute charge collector. In addition, calculations have been performed to determine the 
magnitude of several possible sources of error in using the cup. 

SYMBOLS 

df. distance through volume defined by electrodes 

dE/dx rate of energy loss per unit path length in absorber 

e number of delta rays per proton produced in Faraday cup base that escape 

cup 

N,j w number of delta rays per proton scattered into Faraday cup from entrance 

window 

W mean energy required to produce an ion pair in helium 

X distance of lead from point of production of gamma ray to surface of Faraday 

cup 

ji mass attenuation coefficient 

p density of helium 

Subscripts: 

exp experimental 

th theoretical 

FARADAY CUP DESIGN AND CONSTRUCTION 
Cup Design 

During the initial planning for the Faraday cup, it became apparent that (1) one cup 
at a fixed location in SREL could not meet the varied requirements of the laboratory users, 
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(2) constructing three or more Faraday cups for the various experimental areas of SREL 
was not practical, and (3) a movable self-contained Faraday cup system would provide the 
greatest efficiency and utilization. Because no known Faraday cup had been constructed 
for 600-MeV-proton charge collection, the author and SREL staff members decided arbi- 
trarily to design the base of the cup with a thickness at least 15 percent greater than the 
range of 600- MeV protons and to have a reentrant cup with the reentrant portion at least 
four times the inside diameter of the cup. The proposed diameter of the cup, 40.6 cm, 
would be two times the beam -tube diameter of the beam transport system. The estimated 
weight of a cup of these dimensions was over 1360 kg including the vacuum system. The 
large weight and size of the proposed cup would make it too difficult to use or to move 
from one area to another. Thus, the decision was made to design the cup for use with 
proton beams of 10-cm diameter or less, since beams of this size were easily obtained 
at that time. 

After a maximum beam diameter was selected, the outer diameter of the cup was 
made equal to twice the beam diameter or 20.32 cm. The inner diameter of the reentrant 
portion of the cup was made 15.24 cm to give a wall thickness of 2.54 cm. The length of 
the reentrant portion was chosen to be five times the inner diameter or 76.2 cm. Since 
the range of 600-MeV protons is 25.2 cm in copper and 25.99 cm in lead, the depth of the 
cup base was made approximately 17 percent greater than either range (30.48 cm). The 
outside diameter and the length of the cup were set at 20.32 cm and 106.7 cm, respectively. 

The first choice for the cup material was copper. However, the time and cost 
required to produce a copper cup were considered excessive. The cup was then designed 
to have a lead core with a copper jacket. Thus the cup alone weighed 235.9 kg. 

After construction of the cup had begun, the decision was made to build an electro- 
magnet to suppress delta rays from the entrance window and electrons back streaming 
from the cup base. A cylindrical thin- walled (0.32-cm-thick) copper extension was 
designed to attach to the cup, and the vacuum housing was modified to accommodate the 
extention. The magnet was designed to produce a 0.1-tesla field across a 15.25-cm air 
gap. The vacuum housing extension was designed to extend beyond the cup extension so 
that the cup extension was slightly more than half way across the magnetic volume. This 
configuration was chosen so that electrons back streaming from the cup would be deflected 
into the wall of the cup extension and the delta rays from the entrance window would be 
deflected into the wall of the vacuum housing. 

Cup Construction 

A schematic diagram of the Faraday cup assembly is shown in figure 1. The shell 
of the cup was constructed of a copper sheet 0.32 cm thick. The primary absorber was 
commercially available lead which was poured into the copper shell. All seams in the 



jacket were silver soldered to hermetically seal the lead and prevent possible outgassing. 
The interior of the jacket was thoroughly cleaned prior to the lead fill, and then all exte- 
rior surfaces were cleaned to remove dirt and oil acquired during construction and 
handling. 

The three insulating supports for the cup were machined from tetrafluoroethylene 
rods, 30.48 cm long and 7.62 cm in diameter, and each support had a saddle with a radius 
of curvature equal to that of the cup. The minimum thickness of each support was 
5.08 cm. Each support was thoroughly cleaned and then bolted to the vacuum chamber 
base plate. (Minimum cup to bolt distance was 5.08 cm.) The cup was then mounted on 
the insulating supports with the three 1. 91- cm- long copper pegs on the cup placed in the 
predrilled hole in each support. The pegs prevented longitudinal slippage of the cup and 
the cradle minimized lateral motion. 

The estimated weight of the Faraday cup system, 616.9 kg, is given in table I along 
with the weights of the major components. This estimate does not include such items as 
bolts, ion gage readout controller, electrical connectors, and support brackets for the 
magnet. Initially there was come concern that the weight of the cup would deform the 
supports and cause the cup extension to contact the vacuum chamber. However, after 
7 yr of usage, there has been no noticeable change in the supports. 

Electrical connection to the cup was made with a bare copper wire, 0.16 cm in diam- 
eter, which was bolted to the cup and soldered to a coaxial feed through in the vacuum 
chamber base plate. 


Vacuum System 

The vacuum chamber for the Faraday cup was constructed from stainless steel with a 
0.00162-mm finish on interior surfaces. All feed throughs were made in the 1.9-cm-thick 
base plate which was bolted to the upper removable portion of the chamber. An O-ring, 
0.35 cm in diameter, was used to hermetically seal the upper chamber to the base plate. 
The chamber extension was constructed of aluminum and bolted to the chamber. The 
entrance window to the cup was 0.05-cm-thick aluminum foil. 

The pumping system consisted of a 750 liter/sec oil diffusion pump equipped with a 
liquid- nitrogen cold trap and a 0.00613 m^/sec mechanical pump. The diffusion pump was 
mounted directly to the chamber base plate and pumped the chamber through a 10.16-cm- 
diameter hole in the plate. With the diffusion pump hot and 298 K water to cool the pump, 
the chamber could be pumped from atmospheric pressure to less than 5 x 10 torr 
(1 torr = 133 Pa) in approximately 10 min. After approximately 1 hr the pressure would 
reach a steady level of 10 to 2 x 10"® torr. When liquid nitrogen was used in the cold 
trap, the chamber pressure could be held at less than 5 x 10'® torr. 
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Two major problems related to ionization of the residual gas were discovered and 
corrected during the initial operation of the cup. In the original design of the Faraday 
cup assembly, it was decided to use a 120 liter/sec titanium sublimation pump, rather 
than an oil diffusion pump, to avoid possible problems with the back streaming of oil. 
However, during the first attempt to use the sublimation pump, large random variations 
in output current were observed which were quickly traced to the back streaming of ions 
from the pump. The magnitude of the variations precluded any use of the cup for accurate 
current measurements. An attempt was made to suppress the ions by installing an insu- 
lated wire grid between the pump and vacuum chamber and by applying a voltage to the 
grid. It was possible to minimize, but not eliminate, the problem. The reduced flow of 
ions still produced unacceptably large current variations. The titanium sublimation pump 
was then replaced by a 400 liter/sec oil diffusion pump equipped with a baffle and cold 
trap. 

The second problem was related to the use of an ionization gage to measure the 
chamber pressure. The location of the ionization gage was directly below the front end 
of the cup, and the cup intercepted the ions streaming from the gage filament. The 
10" ^ torr range setting of the ionization gage control unit produced a steady 1-juA current 
output. With the gage off, random current fluctuations no greater than 0.01 pA were 
observed. 


Faraday Cup Carriage 

An open rectangular box frame (101.6 cm high, 45.72 cm wide, and 137.16 cm long) 
was constructed from 5.08- by 5.08- by 0.32-cm steel angles. All corners were rein- 
forced with 0.32 -cm -thick gusset plates welded to the angles. The frame was then welded 
to two 111.76-cm-long sections of 15.24-cm steel channel. Two 360° swivel caster 
assemblies (steel wheels 12.7 cm in diameter) were mounted to each channel section, and 
leveling jacks were mounted outside the casters on the channels. The base plate of the 
vacuum chamber was bolted to the frame to increase the strength and rigidity of the frame. 
Lifting the assembly by using the leveling jacks at opposite corners produced no apparent 
twisting of the frame. 

Positioning of the cup along the beam center line has proved to be relatively simple 
with the aid of a plum bob and level. After some experience in using the cup, most exper- 
imenters could position the cup, attach cooling water lines to the diffusion pump, make 
the necessary electrical connections, and pump down the vacuum chamber to less than 
5 x 10" 5 torr in approximately 1 hr. 

Magnet Design and Performance 

A C-type electromagnet was designed and constructed for use as a delta ray sup- 
pressor at the mouth of the cup. The magnet yoke was constructed of SAE 1010 steel, and 
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the magnet coil was made from B & S 17-gage copper wire wound on a copper spool. A 
12.70'cm-diameter shaft was fitted into the spool end, and the arms of the magnet were 
bolted to the shaft. The pole faces were tapered from 17.78 cm to 13.97 cm in diameter 
with an air gap of 16.51 cm. Total weight of the magnet was approximately 81.6 kg. 

A regulated 1.5-kW dc power supply was the current source. The field at the cen- 
ter of the air gap, as measured with a Hall probe, was 0.11 tesla for a current of 3 A at 
170 V dc. Figure 2 shows the cup extension, and figure 3 shows the magnetic field 
strength along the extension center line for 3- A current in the coil. The reduction in the 
maximum field strength was attributed to shielding by the 0.64-cm-thick aluminum wall 
of the vacuum tank extension and parasitic loss in the nearby steel vacuum tank. The 
shape of the intensity along the center line is related to the magnet face design. 

Because the magnet coil was not constructed with cooling capability, the coil tem- 
perature and, consequently, the coil resistance limited the duration of magnet operation. 
With the magnet coil at room temperature and a constant current of 3 A, the voltage of 
the constant current supply would increase with the increase in coil resistance from 
approximately 170 V dc to the maximum output of 400 V dc in approximately 30 min. At 
this point the current would begin decreasing as the coil temperature continued increasing, 
and the magnetic flux in the cup extension would decrease with the current. 

As shown subsequently, the delta rays from the entrance window are estimated to be 
the largest source of error. Although this error is calculated to be less than 1 percent, 
an attempt was made to see whether the magnetic field across the cup extension would 
suppress this source of error. The test involved measuring the current from a helium- 
filled ion chamber and from the Faraday cup with and without power to the electromagnet. 
Approximately 25 different beam currents were used for each test. The ratio of ion 
chamber current to Faraday cup current gave the number of ion pairs per proton for each 
current setting. The average value and standard deviation of the ratio were 34.93 ± 1.73 ion 
pairs per proton with 3-A coil current (approximately 0.07 tesla at the center of the cup 
extension). Without the magnet current the ratio was 36.3 1± 1.08 ion pairs per proton, an 
increase of 4 percent from the value with the magnet current. Although the difference 
between the ratios indicates that the magnetic field does suppress the delta rays, the large 
standard deviations for the ratios indicate that this difference may not be statistically sig- 
nificant. In addition, these data were taken on different days with different helium supplies 
and, therefore, the difference between ratios may also indicate a systematic error. 

CALCULATED CURRENT MEASUREMENT ERRORS 
Delta Ray Escape 

The charge measured from a Faraday cup is always in error to a degree because 
of interactions of the incident particle beam with the nuclei and atomic electrons in the 
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entrance window and in the cup base. These interactions generate both fast and slow 
electrons which modify the true incident charge. Delta rays (or knock-on electrons) are 
produced through collisions of the incident particles with atomic electrons. Delta rays 
produced in the entrance window are predominantly scattered forward, and a fraction of 
these electrons enters the cup and reduces the true charge. In appendix A the number of 
delta rays produced in the entrance window and the number produced in the cup base are 
calculated. Because of the reentrant design of the cup, the number e of delta rays 
produced in the cup base that can backscatter out of the cup is much less than the number 
Nd )W °f delta rays scattered into the cup from the entrance window. The calculated per- 
cent decrease from the true proton current is 


■^d,w “ ^d,e 
100 


0.88 percent 


Tertiary Electron Production 

Tertiary electrons are low-energy electrons produced in electron- electron colli- 
sions. Relativistic electrons from the entrance window and cup base interact with the 
various parts of the cup assembly to give rise to this additional source of electrons. 

Since this is a tertiary process, the number of electrons produced is much smaller than 
the number of delta rays. 

Since there was no electrical insulation between the diffusion pump and the vacuum 
chamber, the chamber and entrance window were grounded through the diffusion pump; in 
addition, a second ground was used to insure a zero potential on the vacuum chamber. The 
electrometer used to measure the cup current was of the vibrating reed type and main- 
tained the cup potential at a few microvolts above ground potential. This slight positive 
potential on the cup not only tended to attract low-energy electrons from around the out- 
side of the cup but also tended to retain the tertiary electrons produced in the cup. The 
net result was a negligible decrease in the measured current from the true current. 

Compton Electron Production 

Another source of error in the measured current from the Faraday cup is loss of 
charge from the outer surface of the cup by electrons ejected from the surface through the 
Compton effect. To estimate the charge lost by Compton electrons, the following assump- 
tions were made: 

(1) Each incident proton undergoes an inelastic collision and produces one 3-MeV 
gamma ray. 

(2) The probability of producing a gamma ray at a given distance in the base is a 
function of the inelastic proton- nucleus cross section. 
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(3) The gamma rays are emitted isotropically from the point of the inelastic 
collision. 

(4) The probability that a 3-MeV gamma ray will reach the copper jacket on the cup 
is given by exp(-jutX) where p is the mass attenuation coefficient and X is the dis- 
tance in g/cm2 G f lead from the point of production to the surface. 

A Monte Carlo computer program incorporating these assumptions was written to 
describe the transport of protons into the cup base, the emission and attenuation of gamma 
rays, and the production of Compton electrons. The results of the Monte Carlo calcula- 
tions are given in table II for different incident beam radii, and the details of the calcula- 
tions are found in appendix B. 


Ion Formation 

The residual gas in and around the cup is ionized by the incident proton beam and, 
to a lesser extent, by secondary electrons and gamma rays. If the gas pressure is suf- 
ficiently high, the resulting ions can cause fluctuations in the observed current; or if 
there is a bias on the cup, there may be a net migration of ions to the cup which would 
bias the current output. As shown in appendix C the number of ion pairs produced per 
incident proton at the normal operating pressure of 5 x 10"® torr is negligibly small com- 
pared with other possible errors. 

Inelastic Proton-Nucleus Collisions 

In appendix D the loss of charge due to protons escaping from the cup through vari- 
ous nucleon-nucleon and proton- nucleus interactions is calculated. The most likely mech- 
anism for proton leakage from the cup is neutron-proton collisions near the surface of the 
cup. In appendix D the details of the calculations are given along with the various assump- 
tions that were made. Table III shows the percent charge lost by escaping protons for 
incident parallel proton beams with diameters from 0 to 12 cm. 

The cross sections for the production of pi-mesons are generally 1 to 2 orders of 
magnitude below the cross sections for nucleon emission around 600 MeV. Therefore, 
the expected charge loss from the cup by mesons is negligible compared with charge 
losses due to neutron-proton interactions. 

Strain Currents 

Strain currents from the tetrafluoroethylene insulators were found to be less than 
0.01 pA which was the lower limit of the electrometer used for the measurements. Gen- 
erally, the Faraday cup is used for beam current measurements above 1 pA and the strain 
currents are ignored. 
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ESTIMATE OF TOTAL ERROR 


Table IV summarizes the estimated error from various phenomena associated with 
the cup. The values indicated are taken for a 12 -cm -diameter input beam as the worst 
case. As shown in the table, the total estimated error was -0.83 percent, and delta rays 
from the window were the largest source of error. The other errors were an order of 
magnitude less than the delta ray error. 

CALIBRATION OF ION CHAMBERS WITH THE FARADAY CUP 

The primary use of the Faraday cup at SREL has been to calibrate ionization cham- 
bers and to measure the full intensity of the external proton and alpha- particle beams. In 
a series of experiments to obtain the differential scattering cross sections for protons of 
various elemental targets, a helium-filled ion chamber (ref. 7) was calibrated with the 
Faraday cup and then used as the beam current monitor during the measurements. The 
ionization chamber was constructed from drawings supplied by the Neutron Physics Div- 
ision of the Oak Ridge National Laboratory. 

The calibration curve for the ion chamber is shown in figure 4. The number of ion 
pairs per proton is given by the ratio 

Ion pairs _ Integrated ion chamber current 
Proton Integrated Faraday cup current 

The instrumentation used to obtain the data was calibrated against a secondary current 
standard which had a stated worst-case accuracy of ±1.6 percent from 0.001 nA to 0.01 nA, 
±1.3 percent from 0.01 nA to 0.1 nA, ±1.1 percent from 0.1 nA to 1 nA, and ±0.8 percent 
from 1 nA to 10 nA. 

The data of figure 4 gave the number of ion pairs per proton as 36.31 with a standard 
deviation of 1.08. The chamber was operated at 100 V, and helium flow through the cham- 
ber was approximately 236 cm3/sec. The helium used was of standard commercial grade. 

The response of the ion chamber to 600- MeV protons can be computed from 

= <^)V = 37 05 

V Proton / W 

' 'th 


where 

dE/dx = 2.4436 M £ V 7 CI H 2 ( re f. 8) 
g- proton 
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W = 29.9 — — 

Ion pair 


dt = 2.54 cm 


p = 0.17847 -HifL 
cm* 5 

The value for W used in the calculation was taken from reference 9 where the number 
was derived from experiments with- 340- MeV protons. No estimate of the uncertainty in 
the value of W was given. A value of W of 30,2 eV is reported in table 6 on page 233 
of reference 10 for polonium alpha particles. 

The experimental value 


Ion pairs 
Proton 


exp 


36.31 ± 1.08 


obtained from the data in figure 4 does not include the calculated systematic error of the 
cup. Using the error bias of -0.83 percent yields a new value of 


/ ion palrs \ =36 . 01±1 . 08 
\ Proton / exp 

which agrees with the calculated value within the uncertainty of the measurement. If this 
value is used to calculate the energy required to produce an ion pair in helium, the result 
is 


W 


exp 


(dE/dx)d t P 

'"ion pair\ 

V Proton / 


= 30.76 ± 0.95 eV 


for nominal 600- MeV protons. 


CONC LUDING RE MARKS 

The utility of a mobile self-contained Faraday cup system for beam current meas- 
urements of 6 00 -MeV protons has been demonstrated during the past 7 yr at the NASA 
Space Radiation Effects Laboratory. The Faraday cup system, which has an estimated 
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weight of 616.9 kg, is mounted on casters and can be easily moved. Normally one person 
can set up the system for beam current measurements in approximately 1 hr. 

Several possible errors in using the cup for beam current measurements of 600- Me V 
protons have been evaluated. The calculated systematic error for all sources considered 
is approximately -0.83 percent. The largest single source of error arises from delta rays 
entering the cup from proton-electron interactions in the entrance window. For nondi- 
verging, incident proton beams, the errors associated with all other radiation -producing 
mechanisms in the cup are an order of magnitude less than the delta ray error. 

From data obtained in calibrating helium-filled ion chambers with the Faraday cup, 
the mean energy required to produce one ion pair in helium was found to be 30.76 ±0.95 eV 
for nominal 600- MeV protons. This value agrees well, within experimental error, with 
reported values of 29.9 eV obtained with 340-MeV protons and 30.2 eV obtained with polo- 
nium alpha particles. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 5, 1975. 
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APPENDIX A 


ERROR PRODUCED BY DELTA RAYS 


To estimate the effect of knock-on electrons from the window, the cross section for 
proton- electron collisions as given in reference 11 is used. This expression, in the nota- 
tion of the present paper, is 


N a Z o 2rn e dT’ 

da = it -2- r „ 2 — £ -£i_ 

A “ 


0 2 (T f ) 2 


1 - 0 < 


\ 2 


T ' * \ 2 

1 max \T + m p 


(Al) 


where 

r e classical radius of electron, 2.8178 fm 

m e rest mass of electron, 0.511 MeV 

/3 speed of incident proton in terms of velocity of light, 0.792 

T r kinetic energy lost by proton in collision 

T kinetic energy of proton before collision 

m p rest mass of incident proton, 938.256 MeV 

T niax maximum kinetic energy loss allowed by conservative laws 

Avogadro’s number, 6.0225 x 10 22 atoms/mol 

Z atomic charge of element, 13 for aluminum 

A mass number of element, 26.98 for aluminum 

The energy T* (as a function of 6) attained by an electron in collision with a 
600-MeV proton is given by equation 3 on page 14 of reference 11, which, except for nota- 
tion, is 
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APPENDIX A 


T f = 2 m, 


p 2 cos 2 e 

H 


m e + Pn + 


1/2 


-|2 


- pp 2 cos 2 6 


(A2) 


where 

p momentum of proton, 1218.98 MeV/c at 600 MeV 

XT 

c velocity of light 

8 angle of electron deflection relative to proton incident direction, deg 

The maximum energy T’ is given by equation (A2) for 6 = 0°. 

max 

As shown in figure 2, the Faraday cup extension has an outside diameter of 13.02 cm 
and is located 11.69 cm from the entrance window. In use of the cup, the incident proton 
beam was parallel and had a diameter of less than 2.54 cm. However, if the diameter of 
the parallel beam is assumed to be 5.08 cm, it is then possible for delta rays in the angu- 
lar range from Q = 0° to 9 = 38,2° to enter the cup. The corresponding energy of 
these electrons would then vary from T^ iax = 1.72 MeV for 8 = 0° to T 1 = 0.67 MeV 
for 8 = 38.2° as given by equation (A2). 

To obtain , the number of delta rays per proton from the 0.0 5- cm- thick window, 

equation (Al) is integrated over the range of T f and then multiplied by N e , the number 
of electrons per cm 2 in aluminum, and by p = 2.72 g/cm^, the density of aluminum, to 
give 



(A3 a) 


N A Z 

N d , w = P N etf - j- 



2m e 

-1 

0 2 

T f 


I3 2 In (T 1 ) t 1 T 1 
T max 2 (T + m p ) 2 


T max = ^ MeV 
T f = 0.67 MeV 


(A3b) 


N d,w = 8 - 86 X 10 ' 3 


(A3c) 
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where 


N e = N A |pt w 

and where the thickness of the entrance window t w is 0.051 cm. 

As stated previously, the lead core of the cup is surrounded by a 0.32-cm-thick 
copper jacket. The incident protons produce delta rays in this jacket that backs catter 
into the reentrant portion of the cup. Now the range of 1.72-MeV electrons in copper is 
given in reference 12 as 1.13 g/cm 2 . If proton -electron collisions are assumed to occur 
at an average depth of 0.557 g/cm 2 , then only those electrons backscattered with energy 
greater than 0.9 MeV can escape from the base. To estimate the number of delta 

rays which occur in the jacket, the right side of equation (A3b) is evaluated from 
T 1 = 0.9 MeV to T* = 1.72 MeV where for copper Z - 29, A = 63.54, and 

IXlaX 

p = 8.92 g/crn^; the result is 
N d b = 9.24 x 10" 2 

If the delta rays are scattered isotropically in 4?r geometry, then the number of back- 
scattered delta rays per proton that can escape the cup is given by 





= 8.3 x 10' 5 


where the inside diameter of the cup d is 15.24 cm and the distance from the base to 
the near end of the cup extension k is 76.20 cm. 

The loss of electrons from the cup would cause an apparent increase from the true 
proton current, whereas the addition of electrons (delta rays from window) would cause 
an apparent decrease from the true proton current. From these calculations the number 
of delta rays w from the window is much larger than the backscattered delta rays 

that can escape from the cup; thus the decrease from the true proton current would 
be approximately 


^d,w ^d,e no 
— 2 2— = 0.88 percent 

100 
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ERROR PRODUCED BY COMPTON ELECTRONS 

In Compton scattering, gamma rays or X-rays undergo elastic collisions with free 
electrons (assumed to be initially at rest). The kinetic energy T e of a recoil electron 
is given by 


T~ = T 


2a 


1 + 2a + (1 + a)^ tan^ <t> 


(Bl) 


where 


energy of incident photon, MeV 


Of 


= T y/"e 


m. 


rest mass of electron, MeV 


<t> 


angle of recoil electron to direction of incident proton, deg 


The maximum energy of the recoil electron occurs when <p = 0, that is, 


( T e)max = T r(l^) 


(B2) 


The average energy per Compton electron is given by 


e° i 


( T *)av = T r¥ 


(B3) 


where e c^ is the average absorption cross section and e cr is the average collision 
cross section. On page 688 of reference 13, cr is given as 

C 3. 


e CT a = 2,rr e‘ 


2(1 + a)^ _ 1 + 3a _ (1 + q)(2q 2 - 2a - 1) 

a 2 (l + 2a) (1 + 2a)2 a 2 (l + 2a )2 


4a ^ 


r. 

1 + a 


3(1 + 2a) 3 \ a 3 


^ + i ln(1+2a) 


cm 2 / electron 


(B4) 
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and e cr is given on page 684 of reference 13 as 


e CT = 2m e 


/ 1 + a 

2(1 + a) 

V a 2 

_ 1 + 2a 


a 


+ — In (1 + 2a) 
2a 


— l + j£L A cmVelectron 
(1 + 2a) 2 


(B5) 


where r e is the classical radius of the electron (r Q in ref. 13). 

To estimate the charge lost by Compton electrons escaping from the outer surface 
of the cup, the following assumptions were made: 

(1) Each incident proton undergoes an inelastic collision and produces one 3-MeV 
gamma ray. 

(2) The probability of producing a gamma ray at a given distance in the base is a 
function of the inelastic proton-nucleus cross section, 

(3) The gamma rays are emitted isotropically from the point of inelastic collision. 

(4) The probability that a 3-MeV gamma ray will reach the copper jacket on the cup 
is given by exp(-juX) where fi is the mass attenuation coefficient and X is the dis- 
tance in g/cm^ of lead from the point of production to the surface. 

The average energy per Compton electron produced by 3-MeV gamma rays is 
1,73 MeV as calculated from equation (B3). Electrons of this energy have a range in 
copper of approximately 1.16 g/cm^, and since the thickness of the copper jacket is 
2.832 g/cm^, only those electrons produced near the surface will, on the average, escape. 
Of course, the energy distribution of Compton electrons varies from a maximum energy 
of 2.76 MeV (calculated from eq. (B2)) to a zero minimum energy; therefore, the shape of 
the energy distribution and the location of Compton electron production determine the 
number of escape particles. The differential cross section for giving an electron a 
recoil energy between T and T + dT is given (in notation of present paper) by equa- 
tion (41) on page 54 of reference 14 as 
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and the number of Compton electrons produced with energy between (Tgj^^ anc * ® in 
an element of thickness Ax^ in the jacket is given by 

/ A Z p( Te ^max d a 

(N e C ). = N AI PAx i J () _dT 


/ \ ^r 


T + 4o? + 1 j + 1 
2 mJ 2« 2 ) a 2 


-3T y 2 - Ty In (T y - Tj 


> (B7) 


2T. 


/ \2 _ / \ . i 

r / „\i 

(Ty - Tj + 2(Ty - Tj - Ty In 

T r( T y - T ) 


2 

+ a 


T + T y In (T r - 


(Te) 


max 


where 

Avogadro’s number, 6.0225 x 10 2 ^ atoms/mol 
Z atomic charge of element, 29 for copper 

A mass number of element, 63.54 for copper 

p density of element, 8.96 g/cm^ for copper 

The number (^g 0 ). of electrons that can escape the jacket is a function of the depth x 

in the jacket at which they are produced. At a given depth there is a minimum electron 
energy (Te( x i))min below which the electrons cannot escape. Thus, the number of 

Compton electrons N e c escaping from the surface for each 3-MeV gamma ray inter- 
acting within the copper jacket can be estimated by the following expression: 


n 6 c = y «) c (Te)max ^ dT 
~ t le/i W, dT 


min 


(B8) 
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where 


y Xj = 2.832 g/cm^ 
i 

If Np incident protons each produce one 3-MeV gamma ray and if the probability of a 
3-MeV gamma ray reaching the copper jacket is exp(-pX), then the total number N e 
of Compton electrons leaving the surface for Np incident protons is given by 

Np 

N e = N e c l exp (-pXjj - (B9) 

j=l 

To evaluate this equation, a computer program was written which used the Monte 
Carlo technique. The base of the cup (345.6 g/cm^ thick) was divided into 300 segments, 
each 1 g/cm2 thick, to approximate the range of 600-MeV protons in lead, 300.7 g/cm^. 

A beam of 600-MeV protons, assumed to be parallel, is incident on the cup along the cup 
center line. Each proton is multiply scattered through intervening segments until it 
comes to a randomly selected slab. At the selected segment the lateral displacement 
of the proton from its incident trajectory is calculated. Next, a random direction is 
selected, the straight-line thickness X to the surface is calculated, and the coordinates 
of the point of intersection with the surface are determined. The probability that a 
proton-nucleus interaction will occur after traversing a distance H in lead from the 
point of entry to the selected segment is given by 

^inelastic 1 exp^-N^ ^ — dxj 


where a is the inelastic proton-nucleus cross section. The probability value is used to 
form a weighted distribution of gamma rays produced in the cup. 

As an aid in plotting the intensity distribution of electrons leaving the surface, the 
cylindrical surface was divided into 100 bands, each 1.016 cm wide, and the base end was 
divided into 10 concentric rings, each 1.016 cm wide. The coordinates of each electron 
leaving the surface were tested to determine which band was used to catalog the particle. 
Note that the cylindrical bands have equal areas, A C yj = tfRcyl w = 329.48 cm^, but the 
concentric rings on the base have unequal areas given by 


A b = tt 


(n + 1)‘ 


-n2 


w 


(n = 0, 1, 2, . . ., 9) 


where w is the width of the bands and R C yl is the radius of the cylinder. 
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For a total of Np incident protons, the number N e of electrons leaving the sur- 
face was obtained from the program by evaluating equation (B9). However, in figure 5 
the number of electrons leaving the surface in each band is plotted as a function of dis- 
tance (band number) along the surface. The four curves shown in the figure correspond 
to incident beam radii of 0, 2, 4, and 6 cm. Each curve is based on 11 000 incident pro- 
ton histories in the computer program. The percent charge loss from the cup by Compton 
electrons is given in table II. 
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IONIZATION OF RESIDUAL GAS 


To estimate the number of ions pairs formed per incident 600- MeV proton, the 
residual gas is assumed to be principally nitrogen at a pressure of 5 x 10" 5 torr 
(1 torr = 133 Pa). The residual gas occupies a volume whose length l from the 
entrance window to the cup is 114.83 cm. The number of ion pairs per proton is given 

by 

Ion pair _ (dE/dx)ZpR 
Proton W 

where 

2 

dE/dx rate of energy loss per unit path length in nitrogen, 2.265 ^£^~ cn V (ref. 8) 

§ 

p density of nitrogen under standard conditions, 1.2504 mg/cm 3 


W energy required to produce one ion pair in nitrogen, 36.3 eV (ref. 10) 


R 


ratio of chamber pressure to standard atmospheric pressure, 


5 x 10~5 torr 
7.60 x 10 2 torr 


Evaluating this expression shows that 


Ion pairs = 5 9 x 1Q - 4 
Proton 

which is negligibly small compared with other possible errors. 
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CHARGE LOST THROUGH PROTON-NUCLEUS COLLISIONS 


The probability of a proton- nucleus collision is quite large for a 600- Me V incident 
proton coming to rest in lead. It is assumed that the probability of inelastic collisions is 
given by the expression 



or, by changing the variable of integration from path length to energy, 


inelastic 1 ‘ exp /' N A ^ A i dE 

\ P d xj 

where 

a total inelastic cross section, cm^ 

A mass number, g/mol 

i particle path length in absorber, g/cm^ 

Avogadro’s number, atoms/mol 
p density of absorber, g/cm^ 

dE/dx rate of energy loss per unit path length in absorber 
Ej incident kinetic energy 


(Dl) 


Evaluating equation (Dl) for 600 -MeV protons stopping in lead shows that, on the average, 
approximately 0.758 percent of the incident particles undergo an inelastic collision. 

If the energy of nucleons within the nucleus is ignored and it is assumed that an 
incident proton interacts with individual nucleons or small clusters (deuterons, tritons, 
helium nuclei, etc.), then from relativistic kinematics, the maximum energy of the most 
prevalent charged secondaries is given by 
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T = 


2m t cos* 2 9 
— - cos^ 6 

f 


02 ) 


For B = 0, 


2m t 



and 



(D3) 


where 



Pj momentum of incident particle, MeV/c 

Tj kinetic energy on incident particle, MeV 

T kinetic energy of recoil particle, MeV 

2 2 

Ej total energy of incident particle, Tj + irtj , MeV 

m-j rest mass of incident particle, MeV 

rrv rest mass of target particle, MeV 

0 . angle of recoil relative to incident particle direction, deg 

Equation (D3) was evaluated and the range in lead of the five most prevalent sec- 
ondaries was calculated by using equation 14 of reference 15. This range was compared 
with the cup base thickness. Table V gives the results of these calculations. Since the 
thickness of the cup base is 345.6 g/cm^, the loss of charge through the cup end by recoil 
particles is remote. 
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The charge lost from the cup by neutron- proton collisions near the surface was 
considered to be the most probable mechanism for the emission of protons from the sur- 
face. For each proton-nucleus interaction in lead there are a large number of neutrons 
emitted. In figure 6 the average number of neutrons per interaction in lead is plotted as 
a function of incident proton energy. The data points were obtained from calculations 
given in reference 16. Most of the neutrons from the inelastic events are produced 
through nuclear deexcitation with energies generally less than 5 MeV. The cascade neu- 
trons can have energies up to the energy of the incident proton and they can undergo a 
neutron-nucleus interaction with the emission of a proton which may have sufficient 
energy to leave the cup. 

To estimate the number of protons lost from the cup, the following method was 
employed in a Monte Carlo computer program. Since the range of 600- MeV protons in 
lead is 300.7 g/cm 2 , the cup base was divided into 300 slabs, each 1 g/cm 2 thick; it was 
assumed that each incident proton lost kinetic energy ATj in the ith slab and was 
deflected through an average angle (9 ) . by multiple scatterings. The probability P 
that a proton-nucleus interaction will occur in the ith slab is given by 


P = 



a dT 
A ldT 
P dx 


(D4) 


where T^ n and are the kinetic energy of the proton entering the slab and the 

kinetic energy of the outgoing particle, respectively; AT^ = T^ - T 

To simplify the computational details, it is assumed that each proton-nucleus colli- 
sion can be replaced by an elastic collision between the incident proton and a neutron or 

cluster of neutrons. The momentum p’ of the scattered proton and the momentum p ! 

p n 

of the recoil neutron or cluster of neutrons are calculated from the following relativistic 
two-body kinematical equations: 


0Ep C cos £p ± 

t i 

(v ) 2 

. 

- y 2 m p 2 (l - f cos 2 e p ) 

l/Z 

P P ” 

yfi - 

)3 2 cos 2 0pj 




(D5) 


(D6) 


23 



APPENDIX D 


and 


where 

m p 

m n 


E, 


El 


7 n 


and 


sin 

tan 6 n rrr 

/3E 

cos 0 D — — 
p p 

rest mass of incident proton 

mass of neutron or neutron cluster 

momentum of incident proton 

total energy of incident proton 

total energy of scattered proton 

angle of scatter of proton 

angle of recoil of neutron or neutron cluster 


E — Ep + m n 


0 - 



Ep + m n 


y = 


1 


1 - j3 2 





E c = 


( m p + m n V 


+ 2T^m n 


1/2 


(D7) 
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The kinetic energy Tp of the scattered proton is given by 




m p 


and the kinetic energy of the recoil neutron or cluster by 





m n 


A computer program, which incorporated equations (D4) to (D7), was written to sim- 
ulate the energy loss and scattering of the incident protons, proton-neutron elastic colli- 
sions, and subsequent neutron-proton scattering. It is assumed that the incident beam of 
protons has a uniform density distribution and has no divergence. Each proton is multi- 
ply scattered through intervening slabs to a randomly selected ith slab. An angle of 
scatter for the proton in the ith slab is randomly selected over the interval 0 < S ^/2. 
Equations (D5) to (D7) are evaluated by using 6p and the total energy Ep of the proton 
in the ith slab to obtain the scattered proton energy Ep and energy of the recoil neutron 
Ej V Equation (D4) is evaluated to give the probability of an inelastic event. The distance 
from the point of collision to the surface is calculated from the scattered proton direction. 
If the range of the scattered proton of energy Ep is greater than the distance to the sur- 
face, the proton is assumed to have escaped. For a total of Np incident protons, the 
number escaping through this mechanism is denoted by Np. 

To account for the multiple production of neutrons, the data in figure 6 for cascade 
neutrons was used, as follows: for protons in the energy range from 600 MeV to 400 MeV, 
assume three cascade neutrons per interaction; from 400 MeV to 50 MeV, assume two 
neutrons per interaction; and below 50 MeV, assume one neutron per interaction. The 
available recoil energy from the elastic scattering of a cluster of two or three neutrons 
is divided among the particles at random. It is further assumed that the recoil neutrons 
travel parallel to the direction defined by 0 n . Each neutron of the initial cluster has a 
probability P" of having a collision with a nucleus at a randomly chosen point along its 
trajectory. The probability P" of such an event is given by 

P" = 1 - exp ^-N a £ £ dxj (D8) 

where the cross section a is understood to be for a neutron-nucleus collision in lead. 

The inelastic proton-nucleus and neutron-nucleus cross sections are nearly the same as 
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the cross sections shown in figure 4 of reference 17. For each neutron-nucleus event, 
the distance to the surface along the incident neutron direction is compared with the range 
of the tertiary proton. If the range is greater than the distance, the proton escapes from 
the cup. The number of protons escaping through neutron- proton interactions is denoted 
by Np. The total number of protons escaping from the cup for N p incident protons is 
the sum of N p and N”. 

The foregoing procedure indicated that an extremely small number of protons were 
lost from the cup. The program was rerun with the requirement that each secondary 
neutron have a knock-on collision with a proton, that is, P ,T = 1. This overestimate is 
used as the percent charge loss shown in table III. 

The foregoing discussion has assumed a parallel incident beam. No calculations 
have been made for a divergent beam. However, if the incident beam is focused to a point 
on the front face of the base and the beam is entirely within the diameter of the cup exten- 
sion, then the angle of divergence is approximately 3°. If the focus point is moved along 
the center line toward the entrance window, no charge is lost because of penetration by 
the incident beam until a point approximately 50.8 cm from the entrance window is 
reached (divergence of 7.1°). At this point some of the incident protons will pass through 
the sides of the cup. 
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TABLE I.- ESTIMATED WEIGHT OF FARADAY CUP 


Component 

Weight, kg 

Cup 

235.9 

Vacuum housing 

131.5 

Carriage 

81.6 

Vacuum diffusion pump 

49.9 

Magnet 

81.6 

Power supply for magnet 

36.3 

Total weight 

616.9 


TABLE II.- PERCENT CHARGE LOSS THROUGH COMPTON EFFECT 
FOR SEVERAL INCIDENT PROTON BEAM RADH 


Beam radius, 

Charge loss, 

cm 

percent incident charge 

0 

0.073 

2 

.061 

4 

.057 

6 

.063 


TABLE III.- CHARGE LOSS DUE TO PROTON LEAKAGE FROM CUP BASE 
[incident beam assumed nondivergent] 


Beam diameter, 

Charge loss, 

cm 

percent incident charge 

0 

0.009 

4 

.010 

8 

.014 

12 

.020 






TABLE IV.- SUMMARY OF ESTIMATED ERRORS 


Source of error 

Charge loss, 
percent incident charge 

Delta rays 

-0.88 

Tertiary electrons 

(a) 

Compton electrons from surface 

0.06 

Ion formation in residual gas 

(a) 

Neutron- proton collisions 

-0.02 

Pi-mesons 

(a) 

Total decrease from true current 

-0.83 


Negligible. 


TABLE V.- MAXIMUM RECOIL ENERGY OF CHARGED SECONDARIES 
FROM INCIDENT 600-MeV PROTONS 


Secondary particle 

Maximum energy, MeV 

Range in lead, g/cm^ 

Proton 

600 

300.7 

Deuteron 

547 

178.2 

Triton 

479 

110.9 

Helium- 3 

479 

27.64 

Helium -4 

421 

36.7 
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Figure 2.- Faraday cup extension and location of electromagnet pole face. 
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Figure 3.- Magnetic field along center line of Faraday cup from entrance window 

to a depth of 25.4 cm. 3-A current. 
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Figure 4.- Calibration data for helium-filled ion chamber. 
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Figure 6.- Neutrons produced per proton-nucleus interaction in lead. Data points 
were obtained from calculations given in reference 16. Solid lines are used 
only to indicate trend. 
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